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We study the extremal dynamics emerging in an out-of-equilibrium one-dimensional Jepsen gas of 
(TV + 1) hard-point particles. The particles undergo binary elastic collisions, but move ballistically 
in-between collisions. The gas is initally uniformly distributed in a box [-L, 0] with the "leader" 
(or the rightmost particle) at X = 0, and a random positive velocity, independently drawn from 
a distribution 4>(V), is assigned to each particle. The gas expands freely at subsequent times. We 
compute analytically the distribution of the leader's velocity at time t, and also the mean and 
the variance of the number of collisions that are undergone by the leader up to time t. We show 
that in the thermodynamic limit and at fixed time ( > 1 (the so-called "growing regime"), when 
interactions are strongly manifest, the velocity distribution exhibits universal scaling behavior of 
only three possible varieties, depending on the tail of cj>(V). The associated scaling functions are 
novel and different from the usual extreme- value distributions of uncorrelated random variables. 
In this growing regime the mean and the variance of the number of collisions of the leader up to 
time t increase logarithmically with t, with universal prefactors that are computed exactly. The 
implications of our results in the context of biological evolution modeling are pointed out. 

PACS numbers: 02.50.-r, 05.90.+m, 87.10.+e, 05.40.-a 



I. INTRODUCTION 

The study of the statistics of the maximum or the min- 
imum of a set of random variables, generally referred 
to as the extreme- value statistics (EVS), is important 
in diverse areas including disordered systems such as 
spin-glasses [l| and directed polymers @, H, 0|, turbu- 
lent flows la], sorting and search problems in computer 
science @, @, HJ , fluctuating interfaces [H, EH ITll fl2 . EH 
EH EH > granular matter [l6[ , growin g ne tworks [171] . and 
models of biological evolution [18|, |l9|, l2 0l. | 2l |. The theory 
of EVS is simple and well understood [22|, [23| when the 
TV random variables V\, V2, • ■ ■ , Vn are statistically un- 
correlated and each of them is drawn from the same com- 
mon parent distribution 4>{V). Their maximum V max = 
max(Pi , V2, • ■ ■ , Vn) is a random variable whose probabil- 
ity density P(VN)dV = Prob[V < V max ^ V + dV,N] 
is known to have a scaling form for large N, 



P(V,N) 



1 



h(N) 



Gi 



V - aj(N) 
UN) 



I, II, III. 



(1) 

Here the subscript i refers to the three types of tails of 
the parent distribution 4>(V), namely: 

(I) a tail decaying faster than a power law, such as 
(f>(V) oc V 13 exp (-V s ) when V —> 00, with S > 0; 

(II) a power-law tail, such as 4>(V) oc V~^ as V — > 00, 
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where (3 > 2; 
and 

(III) a bounded distribution, such as (j)(V) oc (V c — V) 13 
when V -> V- , with p > 0. 

The functions di(N)-s and bi(N)-s are non-universal scale 
factors that depend on the details of <ft(V). However, the 
scaling functions Gi(z)-s are universal [13, [Hj], in the 
sense that they are only of three possible varieties, de- 
pending on the three classes I, II, and III, but are other- 
wise independent of the details of <fi(V). These universal 
functions are known, respectively, as (I) Fisher-Tippett- 
Gumbel, (II) Frechet, and (III) Weibull distributions. 

In contrast to the above case of independent 
identically-distributed (i.i.d.) stochastic variables, the 
EVS is much less understood when there are correla- 
tions or interactions between the random variables. In 
the presence of static interactions, exact results for the 
EVS are known only in few cases, such as for fluc- 
tuating (1 + l)-dimensional interfaces in their steady 



states 
problems [2 
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and for a class of directed polymer 
In this paper, we present exact 
asymptotic results for the EVS in an interacting parti- 
cle system, for which the interactions between the ran- 
dom variables are manifest dynamically. Our system is a 
one-dimensional gas consisting of (N + 1) identical hard- 
point particles that undergo binary elastic collisions. At 
these instantaneous collisions the particles thus merely 
exchange their velocities, while in-between collisions they 
move freely. Due to the simplicity of the dynamics, this 
so-called Jepsen gas often admits analytical treatments 
for various externally imposed constraints, and hence has 
a rather rich history, see e.g. 0, [H, [H, H3, IH ES H3] ■ It 
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FIG. 1: A realization of the trajectories for N — 5 particles 
plus the leader. The thick line indicates the trajectory of 
the leader, and the labelled arrows refer to the successive 
values of its velocity. The points where the leader trajectory 
gets modified through collisions are indicated by the big black 
dots. The thin lines correspond to the trajectories of the other 
particles. 



has also proved very useful in the study of a class of out- 
of-equilibrium problems like the evolution of the "adia- 
batic piston" [5l|, [S3, [H, HH , the Jarzynksi theorem 35 1 , 
and a quasispecies biological evolution model pjj H"" 



211. 



Moreover, the Jepsen gas also turns out to have impor- 
tant applications in spin transport processes in the one- 
dimensional nonlinear a- model (36l . [37j . 

We start from an initial condition at t = where 
the Jepsen gas of (N + 1) particles occupies the inter- 
val [— L, 0] on the real x axis. The extreme-right "zero"- 
th particle, that we shall conventionally call hereafter the 
leader, is initially located at Xq = 0. The coordinates Xi, 
i = 1, N of the N particles to the left of the leader are 
uniformly distributed in the interval —L^x< Xq = 0. 
Thus the gas has a uniform initial density uq — N/L. 
The initial velocities Vi, i — 0, N of the leader and of 
the other N particles are independent random variables, 
identically distributed according to the parent distribu- 
tion (ftiV)- F° r simplicity, we restrict ourselves here to 
the case of positive velocities, i.e., such that <j){V) = 
for V < 0. No boundaries affect the dynamics of the 
system, which is simply a free expansion. As the system 
evolves in time, the particles collide elastically and ex- 
change their velocities. For a given initial condition, the 
system up to time t is fully described by the set of trajec- 
tories {Xi + Vit, i = 0, 1, ... , N}. Each of the particles 
travels along such a trajectory till it collides with another 
particle, and each collision changes its trajectory. In par- 
ticular, the velocity of the leader increases whenever it 
collides with a particle with higher velocity coming from 
its left, see Fig. [TJ For a fixed particle number N, it is 
obvious that if one waits for a long enough time, then the 
leader will acquire the largest velocity of the initial set 
{Vi, i = 0, 1, . . . N}. Once this happens, the trajectory of 
the leader remains unchanged for all subsequent times. 



In the context of biological evolution of quasispecies, 
this model was first introduced and studied in Ref. [19f . 
The population ni(t) of the i-th genotype or species in- 
creases exponentially with time, ni(t) = ni(0) exp(I^t) 
where the effective rate of reproduction Vi ^ defines its 
"fitness". The logarithmic variable Xi(t) = ln[rii(t)] = 
Xi{Q) + Vit can then be interpreted as the i-th trajec- 
tory of a Jepsen gas. At t = the rightmost particle 
with Xq — and velocity Vq is the leader/most fitted 
genotype; however, if Vq is not the maximal of the ini- 
tial velocities, then it will be overtaken successively by 
faster/better fitted genotypes. At each of these overtak- 
ing events the velocity of the leader changes instanta- 
neously by a finite amount, i.e., the fitness of the lead- 
ing genotype increases discontinuously. These overtaking 
events represent thus the so-called punctuation events in 
the general context of evolution The important ob- 
servable here is the number of leading genotype changes 
up to time t, i.e., the number of collisions that the leader 
undergoes up to time t. In particular, the total number of 
punctuation events till the emergence of the eventual "ab- 
solute" leader genotype exhibits universal dependencies 
on the system size N, that were investigated numerically 
in Refs. 0, H3| and recently analytically in Ref. [2l| . 

The goal of this paper is to compute analytically two 
physical quantities of principal interest, namely: 

(A) the probability density of the velocity of the leader 
P(V, t, N) at any time t, and for a fixed number (N+l) 3> 
1 of particles, and 

(B) the mean and the variance of the number of collisions 
undergone by the leader (i.e., the number of leader's tra- 
jectory changes) up to a finite time t. 

It is clear that in this system, as long as iV ^ 1 is finite, 
there is a natural time scale t*(N) 3> 1 that denotes the 
time taken by the final trajectory of the leader to emerge. 
For t > t*(N), the leader velocity does not change any- 
more, and as such the leader follows on the trajectory 
corresponding to the maximum of the initial velocities. 
Therefore, there are obviously two temporal regimes sep- 
arated by the crossover time scale t*(N), namely: 

The stationary regime (t > t*(N) » 1): In this 
regime, the leader's trajectory remains unchanged for all 
subsequent times. The probability density of the leader 
velocity P(V, N, t) becomes thus time-independent for 
t > t*(N) 3> 1, and it is given by the probability density 
Piy, N) of the maximum of the initial velocities. Since 
the initial velocities are i.i.d., P(V,N) satisfies the scal- 
ing form in Eq. {T]), where the scaling function Gi(z) 
has one of the three universal forms (Fisher-Tippett- 
Gumbel, Frechet, or Weibull), depending on the tail of 
4>{V). Thus, one obtains this regime by taking t — > oo 
limit, but keeping N 3> 1 fixed. The interactions between 
particles, which are manifest only dynamically, become 
completely irrelevant in this regime for the velocity dis- 
tribution of the leader. 

The growing regime (1 t <§; t*(N)): In this regime 
the interactions play an important role, and thus the ve- 
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locity distribution of the leader P(V,t,N) is nontrivial. 
Since t <C t*(N), the finitness of the system does not 
affect the leader dynamics; as such, one can study this 
regime by fixing t and taking the thermodynamic limit, 
i.e., N — > oo limit (at fixed no = N/L). One can show 
that the leader velocity distribution approaches an N- 
independent form P(V,i). In addition, for i > 1 (still in 
the thermodynamic N — > oo limit), P(V, t) has a scaling 
form 



P(V,t) 



6i(t) 



where the index i = I, II, III refers to the three types of 
tails of the parent distribution 4>{V) as mentioned before. 
The functions aj(t)-s and 6j(t)-s are non-universal scale 
factors that depend on the details of 4>(V), but the scaling 
functions Fj(z)-s are universal and are of one of three 
types I, II, and III, but are otherwise independent of the 
details of <fi(V). Moreover, the scaling functions Fi(z)-s 
are different from the scaling functions Gi(z)-s in Eq. ((T|) 
that characterize the EVS of i.i.d random variables. The 
exact forms of the scaling functions Fi(z)-s are detailed 
in the next section. 

We have also computed analytically the mean and 
the variance of the number of collisions n c (t,N) un- 
dergone by the leader up to time t in a system with 
(N + 1) 3> 1 particles. In the stationary regime t > 
t*(N) ^ 1, it is already known that the statistics of 
n c (t,N) become time-independent and presents univer- 
sal iV-dependence [HjHiJEj]. More precisely, (n c (N)) w 
£i ln(iV) , i = I, II, III for large N, where the prefac- 
tor £j is universal and its value is specific to each of the 
three classes i = I, II, III. This result was first conjec- 
tured in Ref. (l9j and was later on proved analytically in 
Ref. [H|. 

In this paper, we compute (n c (t,N)) in the grow- 
ing regime (1 C t C t*(N)) and show that the cor- 
responding mean number of collisions is independent 
of N for large N, and grows universally with time as 
(n c (t)) w 7jln(noi), i — I, II, III, where n$ is the 
initial density, and the prefactor 7, is universal and de- 
pendent on the type of tail of <j>(V). Similar universal 
results are also derived for the variance of n c (t, N) in the 
growing regime (see Sees. II and IV for details). 

The paper is organized as follows. First, for ready 
reference, we present in Sec. II a summary of our main 
results. In Sec. IIII1 we compute the velocity distribu- 
tion of the leader and provide exact asymptotic results. 
Section IIVI discusses the statistics of the number of col- 
lusions undergone by the leader up to time t. In Sec. [V] 
we provide, for illustration, the full analytical results for 
a particular choice of <fi(V). Finally, we conclude with a 
summary and outlook in Sec. IVIl To facilitate the read- 
ing, most of the lengthy calculations are relegated to the 
Appendices A-D. 



II. SUMMARY OF THE MAIN RESULTS 

In this Section, we summarize our main results for 

(A) the velocity distribution of the leader, and 

(B) the mean and the variance of the number of collisions 
undergone by the leader up to time t. 



A. The asymptotic velocity distribution of the 
leader 

As mentioned in the Introduction, the asymptotic ve- 
locity distribution of the leader, both in the stationary 
as well as in the growing regime, is universal, in the sense 
that it depends only on the tail of the parent distribu- 
tion <j>(V). As such, three different universality classes 
emerge. We summarize below these universal behaviors 
for these three classes. 

Class I: Tail decaying faster than a power-law as 



^{V) = AV^e- vS 



(3) 



with 5 > and A a constant related to the details (in 
particular, the normalization factor) of <f>(V). In this 
case, the crossover time scales with N as t*(N) oc N for 
large N. 

The stationary regime (t 3> t*(N) 3> 1): the probability 
density P(V, N) of the leader becomes time-independent 
and is given by the EVS of the i.i.d. initial velocities. In 
the limit of large N (see Sec. IVB for details), P(V, N) 
satisfies the scaling form in Eq. |T]) with 



ai (N) 
h(N) 



X/5 



(1-S)/S 



(4) 



and the scaling function is equal to the Fisher-Tippett- 
Gumbel probability density function (p.d.f.), 



Gi{z) = exp [— z — exp(— z)] 



-00 < z < 00. 



(5) 



The growing regime (1 < t < t*(N)): the velocity dis- 
tribution P(V, t) of the leader becomes independent of 
N as N — > 00 (with t fixed), and has the scaling form 
described in Eq. with 



ai{t) 



b^t) 



fn tA\ 
111 [~P~ ) 



In 



notA 



{IS)/ S 



(6) 



and with the universal scaling function 



Fi{z) 



dUt 



-e' z Ei(~e- z ) , 
-00 < z < 00 , (7) 
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FIG. 2: The universal functions Fi(z) (with continuous line), 
Fjj(z) (dashed line), and Ffjj(z) (dashed-dotted line). The 
superscript "4" represents the value of the parameter 0. 



where Ei(z) is the exponential- integral function j38|. The 
profile of Fj(z) is represented in Fig. [2j and its asymp- 
totics are 



(z — C)e z for z — > oo , 
exp(— e~ z ) for z — ► — co , 



(8) 



where C = 0.577215... is Euler's constant. Correspond- 
ingly, the mean asymptotic velocity of the leader in- 
creases with time t as 



<V(t)> « [Hn t)] 1/S , 



(9) 



and naturally this is also independent on other details of 
Class II: Power-law-decaying distributions <f>(V) as 



4>^\V) = BV 



-0 



(10) 



where > 2 in order to ensure the finitness of the first 
moment of </>(V), and B is a constant. In this case the 
crossover time t*(N) fa iV(0- 2 )/(0-i) f or i arge ]y, and the 
leader velocity distribution has the following asymptotic 
universal form in the two regimes: 

Stationary regime: the time-independent P(V, N) satis- 
fies the scaling form in Eq. |T]), with 



an (AO = 0, 



MAO 



BN 

0~T 



1/03-1) 



(11) 



and the scaling function is the Frechet p.d.f.. 



09-1) 



G H (z) 



exp 



1 



,3-1 



for z > , 



for z < . 



(12) 



Growing regime: the leader velocity distribution P(V, t) 
becomes independent of N as N — > co, and satisfies 
Eq. (O with 



an(t)=Q, b H (t) 



not B 



09-1)09-2) 



l/W-2) 



(13) 



The universal scaling function Fu(z) depends only on the 
parameter (but is otherwise independent of the details 
of 000): 

F^(z) = 09 - l)*"" / dU [/-(/3-i)/(/3-2) e -£/ , (14) 



2-/3 



that can also be written as 

F&(*) = (/3- 1)09- 2)*"" 



dt/< 



(15) 



and is defined for 2 S [0, co). For z < 0, Fjj(z) 
trivially. Its limiting behavior is given by 

03 - 1)09- 2) 



*3-l 



for 



CO , 



/9-1 



(16) 



exp(-z 2 ~ /3 ) for z -> 0+ . 



The mean velocity of the leader increases asymptotically 
as a power-law, 



<n*)> 



0-1 
0-2 



B 



09-1)09-2) 



1/(0-2) 



/?- 3 
^2 



x (not) 1 /^- 2 ) 



(17) 



The above expression is finite provided that > 3. Thus, 
although P(V,i) exists for all > 2, its first moment is 
infinite for 2 < /3 < 3. 

Class III: Distribution <^(V") with a finite maxi- 
mum velocity V c , 



i tail 



(V) =C(V C - V f for V < V c , 



(18) 



with > 0. In this case, f(/V) w yy(0+ 2 )/(0+i) for large 
AT. 

Stationary regime: P(V, N) satisfies the scaling form of 
Eq. (HJ with 



a m (N) = V c , 



5/// 



(AO 



C7V 
/3 + 1 



-1/(0+1) 



(19) 



and the scaling function is the Weibull p.d.f., 

for z > . 



G///(z) 



1) M' 3 exp [-N^ 1 ] for z<0. 



(20) 
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Growing regime: In this regime, the leader velocity dis- 
tribution P(V,t) becomes indepedent of N as N — > oo, 
and satisfies to Eq. (jSJ) with 



ani(t)=V c , b IH (t) = 



(/3+l)(/3 + 2) 



1 1/09+2) 



notC 



(21) 

The /3-dependent universal scaling function is given by 



FL{z) = (/3+l)(/3 + 2)|zp 



dUt 



kl 



(22) 

for z G (—oo, 0], and Fj II (z) = for z > 0. Its limiting 
behavior is given by 

f 03 + 1) 



M 

(/3 + i)r 



exp[— |z|^ +2 ] for z^— oo, 
1 



13 + 2 



\zf for 



0" 



The asymptotic mean velocity of the leader approaches 
the maximum allowed value V c as a power law, 



(V(t)) 



« v; 
x r 



/3+1 //3 + 3 x 



/3 + 1 

' f3 + 2 
13 + 3 
/3 + 2 



[p + 1)09 + 2) 



(«o*) 



C 

-1/(0+2) 



1/08+2) 



(24) 



B. Collision statistics of the leader 



We have also computed analytically the mean and the 
variance of n c (t, N), the number of collisions undergone 
by the leader till a given time t and for a given N 3> 1. 
As before, one is led to consider the two regimes, namely 
the stationary regime (t > t*(N)) and the growing regime 
<**(#)). 

In the stationary regime the mean and the variance 
become time-independent, and they increase logarithmi- 
cally with N for large N. The mean behaves as 



(n c (N)} «&ln(AT), i = 1,11,111, 



(25) 



where the universal prefactor £j depends on the three 
classes of 0(V), 



, 1 



/3-2 
2/3 - 3 ' 

/3 + 2 
= 2/3 + 3 ■ 



and 



(26) 



Similarly, the variance in the stationary regime, for large 
N, behaves as 



(nl(N)) 
where 



(n c (JV)) 2 w o-i ln(JV), i = 1,11,111, (27) 



07 



@ ill 



1 

4 " 

(/3 - 2)(2/3 2 -6/3 + 5) 
(2/3-3) 3 

_ (/3 + 2)(2/3 2 + 6/3 + 5) 



and 



(28) 



(2/3 + 3) 3 

The above results were first derived analytically in 
Ref. [2l[, using a different method. This paper provides 
thus an alternative derivation. 

In the growing regime, we show that the mean number 
of collisions increases logarithmically with t, 

(n c (t)) w 7i ln(not), i = 1,11,111. (29) 

The universal prefactor 7$ is characteristic to each of the 
classes of the parent-distribution (p(V), namely 

1 

2 ' 



In 



Km 



(3-1 
2/3-3 

(3- 



and 



(30) 



1 



2/3 + 3 ' 

Moreover, as discussed in Sec. IIV1 one can also infer the 
variance of the number of collisions, 

(nl(t)} - (n c (t)) 2 ~ ViHnot) , % = 1,11,111, (31) 
where 

1 

m = v 



mi 



(/3- l)(2/3 2 -6/3 + 5) 
(2/3-3)3 

_ (/3+l)(2/3 2 + 6/3 + 5) 



and 



(32) 



(2/3 + 3) 3 

The fact that {n 2 c (t)) — (nAt)) 2 ^ (n c (t)) indicates, con- 
trary to previous claims [33j , that the collision process in 
the thermodynamic limit of the Jepsen gas is not Pois- 
sonian. 



III. VELOCITY DISTRIBUTION FUNCTION 
OF THE LEADER 

A. General relations 

In order to establish the characteristics of leader's 
stochastic dynamics, we shall follow the type of reasoning 
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and the convenient notations of Ref. [32| . Recall that the 
particles simply exchange velocities upon collisions and 
cannot move across each other; therefore at any time t 
the leader rides the instantaneous rightmost trajectory 
among the free trajectories {Xi + Vit, i — 0,1, N}. 
As the number of particles on the left-hand side of the 
leader is a conserved quantity, we can identify this instan- 
taneous trajectory (X p + V p t) of the leader by imposing: 

N 

0(X p + V p t-X t -V l t) = N, (33) 

2—0, i^p 

8 being the Heaviside step function. 

These elements are sufficient to find the conditional 
probability distribution of the leader Pl(X, V,t\Vo) at 
time t > as 

N 

P L (X, v, t\v ) = (Y, S(X -X p - Vpt) S(V - V p ) 



p=0 



N 



S Kr ( N, J2 9 ( X p + V P l - X -- V J) ] ) > 

i—0, i^p 



(34) 



where the brackets (...) denote averaging over the initial 
positions and velocities of the gas particles in [— L, Q], 
The initial condition is obviously 

P L (X,V,t = 0\V ) = 6(X)6(V-V ). (35) 

As shown in Appendix A, one finds 

dz 



P L {X,V,t\V Q ) = j 



2niz 



[A(z, X/t | L/t)} 



N 



x | S(X- Vt)S(V - V Q ) 
+ n (f>(V)9(Vt - X)6(X -Vt + L) 

1 + (Z - l)0(Vqt - X) 

X A{z, X/t | L/t) 
Here T is the unit-circle in the complex z-plane, and 



(36) 



A(z, X/t | L/t) = l + (z-l) 



(L+X)/t 



dUcj>{U) 



f (L+X)/t , x 

dU U 



L 



x/t 



4>(U) 



(37) 



with tiq — N/L the initial density. One can also check 
the normalization of Pl(X, V,t\Vo), see Appendix A. 

We now focus on the growing regime, i.e., we fix the 
time t and take the thermodynamic limit L — > oo, N — > 



oo at constant no- One finds then for the conditional 
probability distribution of the leader: 



P(X,V,t\V ) 
dz 

exp 



2iriz 



not(z 



1) / dU 
'x/t 



U 



X 

T 



cj>{U) 



x {8{X - V t) 5(V - V Q ) + n <l>(V)6{Vt - X) 

x [l + (z-l)9(V t-X)}} . (38) 

Introducing the function 

/•OO 

a(W)= dU (U - W) <j>{U) , (39) 
iff 

one obtains finally: 

P(X, V, t\V ) = e - n ° ta ( V ^ S(X - V Q t) 5(V - V ) 
+n o 0(V) e ~ n Gta{X/t) e{ y t _ _ y ot ^ (4Q) 

This result has a simple physical interpretation, based 
on a "flux argument" . Let V be the instantaneous veloc- 
ity of the leader at time t. The leader's trajectory can 
get bypassed only by the trajectories with higher initial 
velocities U > V; then the rate at which the trajectory 
of the leader gets bypassed is proportional to the flux 
of particles trajectories of slopes higher than V. This 
flux is clearly proportional to the particle density no and 
to the relative velocity (U — V)6(U — V). Therefore, 

/•OO 

noaiy) = no dU(U — V)4>{U) represents the total 
Jv 

instantaneous rate at which the leader's trajectory with 
velocity V gets bypassed by other trajectories. 

The first term in the r.h.s. of Eq. ([40]) represents the 
probability that the leader's trajectory has never been 
bypassed till time t. One notices that this probability 
is exponentially-decreasing with time. In this case, the 
final velocity is Vq and the final position is Vot, thus ex- 
plaining the two 5 factors. Indeed, as seen above, in an 
infinitesimal time dt the trajectory of the leader with ve- 
locity Vo gets hit with probability noa(Vo)dt, and so it 
does not get hit with probability 1 — noot{Va)dt. As such, 
the probability for the leader to keep its initial trajectory 

till time i is [1 - noaiV^dt] 1 1 dt -> exp[-n ta(V )]. 

Let us turn now to the second term of the r.h.s. of 
Eq. (|4T))) . which takes into account all the situations when 
the velocity of the leader got modified through collisions. 
The leader is at X at time t, with velocity V. For this 
event to happen, no trajectory must hit the straight line 
of slope X/t till time t, see Fig. [3J All possible trajecto- 
ries of the leader must lie below this line, and (according 
to the argument above) this happens with probability 
exp[— nota(X/t)]. Now, out of all the trajectories satis- 
fying this criterion, we are interested only in those that 
actually hit the line of slope X/t exactly at time t, and 
exactly with velocity V; this fraction is nQ(/)(V). Thus, 
the total probability for the leader to be at X with ve- 
locity V is no4>(V) exp[— nota(X/t)} (the two 9 functions 
are obvious). 
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X 

X 




Here /i > is the first moment of <^(V), 



FIG. 3: Illustration of the "flux argument" in deducing the 
conditional probability distribution function P(X,V,t\Vo) of 
the leader in the thermodynamic limit (see the main text). 



The conditional coordinate distribution for the leader 
is obtained by integrating P(X, V,t\Vo) over V: 

P{X,t\V ) = e _n ° to ( y o) 6(X - V Q t) 

+ noe _noto i X /*) / dV^(V)6{X -V t) 
Jx/t 

= ^{e- n o ta (X/t)e(X-V t)} , (41) 

and a detailed discussion of its long-time properties, cor- 
responding to the diffusive regime for the particle, can 
be found in Ref. [27|. 

We now consider exclusively the stochastic behavior 
of the velocity of the leader. The conditional velocity 
distribution is obtained by integrating P(X, V, t\Vo) over 
X: 

P(V,t\V ) = e~ n o ta ( V o) S(V -V ) 

+ n t<t>{V) / dW e- n ° ta ( W U(V -V ). 

JVq 

(42) 

Moreover, by averaging over the initial velocity Vq of the 
leader, one obtains the velocity distribution: 

P(V,t) = / dV <f>(Vo)P{V,t\V ) 
Jo 

= 4>{V)e- n * ta( y) 

/•V r-V 

+ n t4>{V) / dVo0(Vo) dWe- n o ta ( W ). 

JO JVo 

(43) 

Based on the properties of the function a(W) (see Ap- 
pendix B), one reaches finally an expression that proves 
to be more convenient for further analysis, 



P(V, t) = <f>{V)e- n ^ + n t<t>{V) 



dW, 



-n ta(W) 



fx= I dVV<j){V). (45) 
la 



B. Long-time behavior of P(V,t) 

Let us inspect Eq. (|4*4")l in the limit t 3> (n /i) _1 . The 
first term in the r.h.s. is exponentially decreasing in 
time, so it can be neglected in this long-time limit. In 
the second term, due to the presence of the exponential 
exp[— nota(W)], the main contribution to the integral 
will come from the W sector for which a(W) is small; in 
view of the property (i) in Appendix B, this happens for 
large values of W, so one can write with a good approx- 
imation: 

/>oo 

P(V, t) « n t(j)(V) / dWexp[-n ta as (W)} . (46) 

Here "as" designates the asymptotic, large- W behavior of 
a(W), which, in view of the definition ([39]) , is determined 
by the asymptotic behavior of 4>(V). We are thus led to 
consider the tail of the distribution <j>(V), according to 
the three classes discussed in Sec. |TTJ 
Class I, Eq. {3|): In this case one finds 

a as (W) fa 4 W p+2{1 - S) exp(-iy 5 ) . (47) 

As shown in Appendix C, this leads to the scaling 
form @ with the scaling function Fj in Eq. and the 
scaling parameters ©• 
Class II, Eq. 03): One has 



a as {W) 



B 



09-1)09-2) W?- 



(48) 



Introducing it in Eq. (|46[) and using a simple change of 
variable 



n Q tB 



L09 -1)09 -2) 



-1/03-2) 



V 



(49) 



one obtains the scaling form © with Fjj as in Eq. (fT5|) 
and the scaling parameters (fl3|) . 

Class III, Eq. U8]): Finally, for distributions with finite 
support 4>(V) one has 



a as (W) 



C 



09 + 1)09 + 2) 
The change of variable 



notC 



09+1)09 + 2) 



l/(/3+2) 



(V c -W) = -z, 



(50) 



(51) 



leads to the scaling form @ with Fjjj given by Eq 



(44) and the scaling parameters in (|2"Tj) . 
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C. Mean velocity of the leader 

Using Eq. |@3) and the properties (B3)-(B6) of a(V), 
followed by a double integration by parts, the mean ve- 
locity of the leader can be written as 



(V(t)) 
+n t 



dV VP(V,t) = fie- 710 ^ 



o 

dV 



a(V) - V 



da(V) 
dV 



-n ta(V) 



Or, with a change of variable 

Z = riQta(V) 
the above integral becomes 



(V(t)> = Lie- not ^- 



The quantity 



notfi 



dZ 



( a(V) 



\\da(V)/dV\ 



V 



(52) 



(53) 



a(V) 



V 



(54) 

has to be expressed as 



\da{V)/dV\ 
a function of Z using Eq. (|53j) . 

In the long-time limit of t 3> (no/j,)^ 1 , the first term 
in the r.h.s. of Eq. (f54|) becomes negligible, and the 
dominant contribution to the integral in the second term 
comes only from the large- V sector. Therefore, 



+ V e" z , (55) 



with Z = noto as (F), so the long-time behavior of the 
mean velocity is determined by the tail of the parent 
distribution 4>(V). Therefore, its expression is the same 
for all the distribution (j)(V) belonging to one of the three 
classes described above, and is given, respectively, by the 
Eqs. ©, OH), and ((24j) above. 

Note that these long-time expressions of the mean 
leader velocity can also be obtained directly by using the 
specific scaling forms J2]) of P(V,t). 



IV. COLLISION STATISTICS 
A. Growing regime 

We are interested in the statistics of the number of 
collisions n c that the leader undergoes till time t in the 
growing regime. As before, we then fix the time t and 
take first the thermodynamic limit N — > oo, L — > oo 
keeping the density n fixed. If V is leader's velocity at 
time t, then the probability that it undergoes one colli- 
sion in the small interval dt is noa(V)dt, see Sec. IHI Al 
Taking the average over the distribution P(V, t) of V, one 
obtains: 



d{n c {t)) 
dt 



n I dV a{V)P(V,t) = n {a{V)) . (56) 



The long-time behavior is obtained by considering in the 
above equation the asymptotic form (0) of P(V,t). One 
notices that the most important contribution to the in- 
tegral over V comes from the region V ~ cii(t),i = 
I, II, III. For the three classes of parent distributions 
4>{V) one obtains, respectively: 
Class I. 



d(n c (t)} 



dt 
n 



S 2 bj(t) ri \ bj(t) 



1 f°° 

-2/ dV e*v[-V s + {13 + 2-28) \nV- 
" Jo 

v-a T (ty 



1 f°° 1 

- / dz e~ z FAz) w — 
tJ-00 2t 



(57) 



Class II. 

d(n c (t)) 
dt 

Class III. 

d(n c (t)) 
dt 



dz z 2 -^F^(z) 



/3 — 1 \ 1 



2/3 -3 I t 



(58) 



1 



/ y ' ' \2(3 + 'd) t 

(59) 



These results lead obviously to the oc ln(not) asymptotic 
growth of (n c (t)), as described by Eqs. (j2"9> J5D|) . 

The mean square number of collisions till time t is given 
by the following integral expression 



P-t pt pOC pOO nV-2 

2 dtx dt 2 / dVi / dV 2 / dU 
Jo Jti Jo Jvi Jvi 



x n (U - V^iU) P(V 2 ,t 2 \U, ti)noa(V 2 ) 



(60) 



The kernel of this five-fold integral corresponds to: 
(i) having one collision in the interval dt\ around t\, pro- 
vided that the leader has velocity V\ at t\ ; (ii) as a result 
of this first collision, an instantaneous change in leader's 
velocity from V\ to U > V\ ; (iii) a second collision in the 
interval dt 2 around t 2 , provided that the leader has veloc- 
ity V 2 ^ U at t 2 (t 2 > ti) and that it had the velocity U 
at ti . The corresponding conditional probability density 
P(V 2 ,t 2 \U, ti) can be easily computed using Eq. (|40|) and 
integrating over the spatial coordinates, and it is found 
to be: 

P{V 2 ,t 2 \U,h) = e " n °( t2 _ i i)«( C/ ) 5(V 2 - U) 
+ n {h-t x )(j>{V 2 ) / dWe~ n o(^ -timW) e( y 2 _ 
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Finally, one has to average over all the realizations of 
leader's trajectory that fulfill conditions (i)-(iii). 

From Eqs. $U§, flST)), and P9"] ) -(|50] ) . one can infer the 
long-time behavior oc ln(noi) of the variance of Tic ; 
resumed in Eqs. (f3Tj) - ([32"]) . This calculation is rather te- 
dious and we have carried it out explicitly only in a spe- 
cial case when (j)(V) — exp(— V) (see Sec. V). However, 
the results for other cases can be inferred by matching the 
late-time growing regime results with that of the station- 
ary regime results derived in Ref. [2l[ at the crossover 
time t — t*(N). This is explained in detail in the next 
subsection. 



B. Stationary regime and the crossover time 

As already mentioned in the Introduction, the Jepsen 
gas for a finite N is completely equivalent to the model 
of an evolutionary dynamics for a quasispecies model in- 
troduced in [HI [2(| • The collision statistics of the leader 
in the stationary regime has been studied before both 
numerically [l^, [2(| and analytically [2l[ , and the results 
in Eqs. (|25|) and (|27|) were derived. Here we will derive 
these results through a different method, that will also 
allow to estimate the crossover time t* (N) for the three 
classes of <j>{V). 

Recall that the crossover time t*(N) is the time at 
which the leader aquires its final, maximum velocity. 
So, at t = t*(N) the typical time-dependent velocity of 
the leader in the growing regime matches with the typ- 
ical value of Vmax = max(Vo, Vi, Vn). To estimate 
,typj we recall that the Vi-s are i.i.d random vari- 
ables each drawn from <j)(V), and therefore, for Af> 1: 



time, 



Prob(V max < V) = f]Prob(Vi < V) 



dU<j>{U) 



i=0 
N+l 



exp 



-N / dU<f>{U) 
Jv 



1- / dU<j>(U) 
Jv 

= exp[Na'(V)} . 



N+l 



(62) 



where the function a{z) is defined in Eq. (|39p . 

Class I: In this case, Eq. (|62p generates the cumulative 
of the Fisher-Tippett-Gumbel p.d.f., 



Prob(V max < V) 



exp 



exp 



V-~a I (N) 



(63) 



where di(N) and bj(N) are given in Eq. ([?]). So in this 
case the typical value of V max is V ma ,x,typ ~ 5/(JV). On 
the other hand, as indicated by Eqs. (T3|), ©, ©, the typ- 
ical leader velocity at time t is roughly Vt yp (i) rj aj(t). 
Matching V maXi ty P = Vt yp (t*(N)) gives the crossover 



t*(N) « —N. 
n 



(64) 



Class II: In this case, one obtains the cumulative of the 
Frechet p.d.f., 



Prob(V max V) w exp 
and hence 



(3-1 



V 



1-/3 



V f w atV(/3-i) 

v max,typ ~ iV 



(65) 



(66) 



On the other hand, according to Eqs. ©, [fig]), ([T7)l . the 
typical leader velocity 



^typ(t) 



Brust 



1 1/03-2) 



.09-1)09-2). 
Then the crossover time is estimated as: 
t*(N) w AT^" 2 )/^- 1 ) . 



(67) 



(68) 



Class III: Finally, one has in this case the cumulative of 
the Weibull p.d.f., 



Prob(V max V] w exp 



JT+T 



(V c - vf +1 



indicating 



V - V , re /v~ 1/(/3+1) 



(69) 



(70) 



The typical leader velocity is given, through Eqs. ((2]), 
(EED, and ([21, by 



v c -v typ (t) 



Cnat 



1 -l/(/3+2) 



.09 + 1)09 + 2) 
and thus the crossover time is 

t*(N)fz iV (/3+2)/(/3+1) 



(71) 



(72) 



Substituting Eqs. JMD, dBSJ, (HD in Eq. one can 
then compute the saturation value of the mean collision 
number (n*(N)) = (n c (t = t*(N))) as a function of N 
as stated in Eq. (f2"5|) . We thus recover, through this 
alternative method, the results derived in Ref. [21] . 

Concerning the variance, the results in Eqs. |27|) . (|28]) 
for the stationary regime were derived exactly in Ref. [2l[ . 
In order that the variance in the growing regime matches 
with that in the stationary regime at the crossover time 
t = t*(N), it follows immediately that in the growing 
regime the variance at late times must behave as in 
Eqs. (|3"Tj) , (f3"2"| . We have verified this prediction by di- 
rect calculation for the special case 4*(y) — exp(— V) (see 
the next Section). However, although desirable, a direct 
calculation of the variance in the growing regime for the 
other cases, using the method outlined in the previous 
subsection, seems too tedious. 
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V. EXACT RESULTS FOR ALL TIME IN A 
SPECIAL CASE 

A particular example that can be studied in full ana- 
lytical detail for all t is 

<P(V)=e~ v (V>0), (73) 

which pertains to Class I with 8=1. By straightforward 
calculations one obtains from Eqs. (|4"4"|) . (f52"| . and (f56")> . 
respectively: 

(i) The probability distribution function for the velocity 
of the leader 

P(Y,t) = e-( n ° t+v ^ <+n te- v [Ei(-n i) - Ei{~n te~ v )] , 

(74) 

with the scaling form ^ corresponding to the func- 
tion |7]) with the parameters 

ai(t) =ln(n t) , 6/ = 1 . (75) 

(ii) The mean velocity is therefore 

(V(t)> = In(no<) + 1 + C- Ei(-not) , (76) 

which for long times is dominated by the logarithmic 
term. 

(iii) The mean value of the number of collisions the leader 
undergoes till time t: 

(n c (t)> - \ [ln(not) +C- Ei(-n t)] , (77) 

with the logarithmic asymptotic increase. 

(iv) Finally, the calculation of the variance of the number 
of collisions is rather lengthy (see Appendix [Djl . but its 
long-time behavior is simply given by 

(n*(f)>-<n c (f)> a «iln("o*), (78) 
in agreement with Eqs. (|3"Tj) and (f3"2")) . 

VI. CONCLUSIONS 

In this paper we have studied the extremal dynamics 
in a one-dimensional Jepsen gas of (N + 1) particles, ini- 
tially confined in a box [-L, 0] with uniform density and 
with each particle having an independently distributed 
initial positive velocity drawn from an arbitrary distri- 
bution 4>(V). We have computed analytically the veloc- 
ity distribution of the leader (or the rightmost particle) 
at time t, and also the mean and the variance of the 
number of collisions undergone by the leader up to time 
t. We have shown that for a given N ^> 1, there is a 
crossover time t*(N) that separates a stationary regime 
(t > t*(N)) from a growing regime (1 -C t < t*(N)). 
While in the stationary regime, the leader velocity be- 
comes time-independent and follows the standard ex- 
tremal laws of i.i.d random variables, it has novel uni- 
versal scaling behavior in the growing regime. The asso- 
ciated scaling functions in the growing regime belong to 



three different universality classes depending only on the 
tail of 4>(V), and they have been computed explicitly in 
Eqs. (O, (fl"5)) . and (J22). These dynamical extremal scal- 
ing functions are manifestly different from the standard 
EVS scaling functions of i.i.d random variables. 

Similarly, we have shown that in the growing regime, 
the mean and the variance of the number of collisions 
of the leader up to time t increases logarithmically with 
t, with universal prefactors that were computed explic- 
itly in Eqs. (|30"|) and (|32f . Also, as a by-product, we 
have provided an alternate derivation of the stationary 
regime results for the collision statistics (mean and the 
variance) of the leader that were obtained in [2l[ through 
a completely different approach. While in this paper we 
were able to compute only the mean and the variance of 
the number of collisions, it would be interesting to com- 
pute the full distribution of the number of collisions up 
to time t, and to compare it with the previously incor- 
rectly suggested Poisson distribution [33j , which remains 
a challenging open problem. 

We have computed here the velocity distribution and 
the statistics of the number of collisions separately in the 
growing regime and in the stationary regime. It would 
be interesting to compute the exact crossover functions 
that interpolate between the two regimes. For example, 
for the velocity distribution of the leader this crossover 
scaling function can, in principle, be computed from our 
general result in Eq. (f3"6")> . which is valid for all times t 
and all values of N. 

The model studied here can also be considered as a sim- 
ple toy model of biological evolution (l9l. I20I l2l| . on which 
some results, numerical and analytical, were known be- 
fore but they were mostly restricted in the stationary 
regime [l9l |20| . [2l| . For example, the number of over- 
taking events of the leader, i.e., the number of punctu- 
ation events till the emergence of the best fitted species 
were studied before [l!| [2TJ, [2l[. However, the authors 
in Ref. [2(3] also studied analytically the distribution of 
the 'label' of the fittest species in the growing regime, 
and by matching the typical leader's label at time t with 
that of the final leader, they were able to extract the 
crossover time t*(N). In this paper, we have studied a 
complimentary quantity in the growing regime, namely 
the distribution of the 'fitness' (velocity) of the fittest 
species. The crossover time t*(N) extracted from both 
of these distributions are in agreement with each other. 
The method presented in this paper may also be useful to 
study the dynamics of other interesting observables in the 
context of biological evolution, such as, for example, the 
persistence of the leader genotype, and the distribution 
of the time interval between two successive punctuation 
events. 
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S(X -X p - V p t)S(V - T/ p ) z -e(*p+^*-*o-v t) 



, FT Z -9(X P + V p t-X 3 -V 3 t)^ 



APPENDIX A: DERIVATION OF EQ. ([36]) AND 
THE NORMALIZATION OF P L (X,V,t\V ) 



We present below the main steps of the derivation of 
Eq. see also [32|]. In order to pursue the calculations, 
it is convenient to use the following integral representa- 
tion of the Kronecker delta: 



= N(5(X -X p - V p t)S{V - V p ))z-^ x - Vo ^ 

x{z -o(x-x J -v J ty- K (A4) 

One has to compute then 



dz 



9(X-X j -V j t)\ 



(Al) 



where p, q are integers, and T is the unit circle centered 
at the origin in the complex z-plane. Using this repre- 
sentation in Eq. (|34| . it results: 



dXj 



dV 3 4>{V 3 ] 



oo 

x-x 



- l + {l-z~ 1 )- j L dX 3 j^dV ]m 



P L (X,V,t\V ) 



x 



z-' + il-z- 1 ) 



dz N 
z 

p ZTTIZ 



(J2 S(X -X p - V p t)S(V - V p ) 

p=0 

N 

x ft z - e ( x p+ v 'p t -^- y J ')} 



dv 3 4>{v 3 : 



dv^iVjWjt-x) 



z^A&X/t | L/t), 



(A5) 



dz 
r 2niz 

N 



X 



S(X - V t)S(V - U )(TT z-HVot-Xj-Vjt)^ where A(z,X/t \ L/t) is given by Eq. (J37J). Note also 
~J[ that 



($2 S(X -X p - V p t)5(V - V p )z- >- 1 1 " 



8{X P +V p t-V t) 



p=l 



N 

x FT z-^Xp+Vpt-x-Vjt)^ 

3=1, j^P 



dX r 



(A2) 



Here &(■■■) is the Heaviside unit step function. Due to the 
initial statistical independence of the particles, one has: 



(S(X -Xp- V p t)5(V - V p )) 

dVp(b(V p )S(X -Xp- V p t)S(V - Vp) 
dX p 8{X -Xp- Vt)cj)(V) 
= y(f>(V)6(Vt - X)9(X — Vt + L), (A6) 
and also 



o 

-L 


-L 



N 



-e^Vot-Xj-Vjt)^ _ u-eiVot-Xj-Vit^N 

= (z-Kx-Xj-WY t (A3) 



z -0(x-v o t) = z -i + (i - z -i)0(y Qt - x) 

= z- 1 [1 + (z - 1) 6{V t - X)} . (A7) 



where the last equality results because of the factor S(X 
V)t) in the corresponding term of the equation (|A2 



Finally, replacing the results (|A3 |) - ([A7)l in Eq. (|A2[) . one 
obtains the expression (|36|) of the conditional probability 
density for leader's coordinate and velocity. One can also 
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check the normalization of this one. Indeed, consider 

/oo roc 
dX / dVP L (X,V,t\X = 0,V ) 
-oo J — oo 

^{ W2 ,V | £ / t )]»4/> 



Using these relations, one can easily obtain Eq. (|44f from 
Eq. m>. 



dV 



APPENDIX C: THE SCALING BEHAVIOR OF 
P(V,t) FOR THE CLASS I OF PARENT 
DISTRIBUTIONS <j>(V) 



tf>(V) [A(z, X/t | L/t)] N ~ l [l + (z- l)6(V t - X)} } 



With the expression (|4T|) of a as , one finds from 
Eq. ([46]): 



at r°° r — t — 
= [A(0,V \ L/t)] N + j- J dX J x dV 

x <P(V)[A(0, X/t | L/t)] N - l 6(X - V t) , 



P(V, t) « ii <^ 



(A8) 



( / dW e 

o 



(CI) 



where for the last equality we used the theorem of 
residues. But: 



d 



dX 
N 



[A(0,V \L/t)} 



N 



[A(0,X/t\L/t)} 



N-l 



dU<t>(U), (A9) 



and therefore: 



AT = [A(0,V I L/t)] N + / dX9(X - V t) 



d 



x^[A(0,X/t\L/t)} N 



[A(0, V | L/t)} N + / dX— [A(0, X/t | L/t)} 



N 



= lim [A(0, X/t | L/t)} N = 1 , (A10) 
which proves the normalization of Pl(X, V,t\Vo). 



APPENDIX B: PROPERTIES OF THE 
FUNCTION a(W) 

Starting from its definition (|39[) . one can easily obtain 
the following properties of the function a(W): 

(i) a(W) is a strictly decreasing function of W 

(Bl) 

(it) a(W) —pt — WforW^O. 

In particular, a(0) = /i, the first moment of 4>(V). 

(B2) 

(in) a(W) -> for W -> oo. (B3) 
(iv) a(W) -» -oo for W -> -oo. (B4) 
da(W) 



(v) 
Note 
(ot) 



dW 
da(W) 



dU<j)(U). 



w 



dW 



w=o 



d 2 a(W) 
dW 2 



<t>(W) 



(B5) 
(B6) 



where 

S(W) = W s - \n(n tA/5 2 ) - (J3 + 2 - 25) In W . (C2) 
Consider 

W = ai (t)+y, (C3) 
where y is small compared to ai(t). Then 
S( ai (t)+y) 

= [ ai (t)] s -ln(n tA/S 2 ) - (/3 + 2 -25) In ai(t) 



(p + 2- 25) 



y + 0(y 2 ). (C4) 



[ ai (t)}i- s aj(t) 
Set 

[a^t)} 6 - ln(n U/(5 2 ) - (/? + 2 - 25) Inaj(i) = , (C5) 
so, to leading order for large t: 

[a/(i)] 5 « ln(n U/,5 2 ) + (/? + 2 - 25)/5 lnln(n U/,5 2 ) , 

(C6) 

(C7) 
(C8) 

(C9) 



and thus 



Then 



where 



ai(t) « [ln(n U/5 2 )] 1/5 . 

W-aj(i) 



S(W = ar(i)+2/) 



bj(t) 



bj(t) « ^(riU/^)]! 1 -*)/ 5 
o 



Now, from Eq. fCT 



P(V,t) w n t6/(*MF^< 



— f 



^g-f^WJ/iutt) 



^/^expI-U 5 +ln(n U) +/31nU] / d[/e 



(CIO) 



where z = (V — ai(t))/bi(t). One has 
expf-U 15 + ln(n U) + (3 In U] 



&?(*) 



e- 2 , (Cll) 



and thus obtains finally Eq. ^ for case I. 
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APPENDIX D: THE MEAN SQUARE 
DEVIATION OF n c (t) FOR cj>(V) = exp(-V) 



where 



Using Eqs. (l60|) for computing (n^t)) For the par- 
ticular case of (f>(V) = exp(— V) (V ^ 0) one has 
a(V) = exp(-U), P{V,t) given by Eq. J74]), and the 
conditional velocity distribution function ()61[) 



P(V 2 ,t 2 \U,h) = e -«o(t2-t 1 )e-^ §( y 2 _ ^ 

+n (< 2 - [JSi(-no(*a - ti)e~ u ) 

-Ei(-n (t 2 - h)e- V2 )}e(V 2 - U) . (Dl) 

Then, after lengthy calculations, Eq. (JSOJ) leads to 

(»c(*)> = ^ ln2 ^+(? + l) 



+ 



7 C 



2C< 



,-n t 



o—n t 



2not 

Ei(-not) — ]n(not)Ei(—not) 
9C C 2 



+ 2e- n ° t Ei(-n t) + (_ + _ + _) + J, 



(D2) 



1 „ f not dr 



'e- noT Ei(-n Q {t - t)) + Ei(-n Q T)] . 

(D3) 



Here V designates the principal part of the above inte- 
gral. 

In the long-time limit n t 3> 1, the main contribution 
to the value of (n 2 c (t)) comes from the logarithmic terms, 



(n 2 c (t))*Un 2 (n t)+(^ + ±) lnl„„M. (D4) 



Combining this result with the expression ([77]) of (n c (t)), 
one obtains finally the asymptotic result (|75|) . 
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